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Effects of perturbations on topologically ordered (TO) states with fractionalized excitations is 
very intriguing. While small perturbations are not expected to affect a robust TO, sizable pertur¬ 
bations must lead to onset of novel routes to order(s). Observability of such emergent orders in real 
condensed matter systems continues to be an open issue of paramount interest. Here, we detail how 
a featureless Kitaev spin liquid with rigorous TO gives way to partial topological ordered state(s) 
coexisting with a wide staircase of novel ‘excitonic’ orders of fractionalized entities in a Zeeman field. 
We provide quasi-exact and novel realization for a range of strange features, from nodal spin-metals, 
fractionalization in selective-Mott states, to Luttinger surfaces and hidden coherence, all character¬ 
istics of the influential resonating- valence-bond as well as fractionalized-Fermi-Liquid views in the 
context of strange metals. As a novel application, we show how Josephson-junction array may 
exhibit fractional Josephson oscillations upon flux tuning as a particularly novel manifestation of 
emergent quantum coherence. 

PACS numbers: 74.70.-b, 74.25.Ha, 76.60.-k, 74.20.Rp 


Topological excitations have fascinated physicists for 
many decades. Remarkably, recent advances in con¬ 
densed matter have reached a stage where the ability 
to observe and even manipulate such exotic excitations 
seem to be within the realm of possibility. Cold-atom set¬ 
ups, topological insulators and appropriately engineered 
Josephson junction arrays (JJA) are poised to offer re¬ 
newed insight into topological excitations arising from 
ground states with topological order (TO). In many in¬ 
stances, however, ensuring their stability is a potential 
hurdle. Another, more exacting hurdle is to unearth 
specific examples of systems capable of exhibiting stable 
(symmetry protected) TO. 

The birth of the exactly soluble Kitaev model (KM) [1] 
led to a spurt of proposals to ‘simulate’ its TO ground 
state in various contexts [2, 3]. In contrast, a micro¬ 
scopic understanding of effects of various perturbations 
is lacking: while mostly studied numerically [4], a co¬ 
herent picture of their effects in terms of changes in 
the spectrum of topological (elementary) excitations of 
the TO state has remained elusive. Given that the KM 
hosts a fractionalized spin liquid phase, novel, unantici¬ 
pated ordered states (partial or complete deconfinement- 
confinement transition(s)), depending upon specific na¬ 
ture of perturbations, are expected theoretically. Classi¬ 
fying such novel orders aids the quest for their observabil¬ 
ity in natural or engineered condensed matter systems, 
apart from the attractive motivation of application(s) to 
fault-tolerant quantum computing. Here, we are moti¬ 
vated by recent proposals in perturbed Kitaev models 
constructed by suitably engineered JJA [3]. Moreover, 
for the realistic range of physical parameters, such JJA 
(see Fig.(l) of supplementary informations) are always 
perturbed by ‘magnetic fields’ along x and/or 2 : bonds. 


Hence, they offer an attractive avenue to investigate 
the above issues in a realistic context. Further, models 
that are variations of the KM and compass models [5] are 
believed to apply even to certain transition-metal oxides 
(TMO) [6], though, in this context, additional perturba- 
tion(s) must be included [7, 8]. In this context, fascinat¬ 
ing issues relate to whether a fractionalized orbital liquid 
can emerge, and whether analogues of orhital-Kondo de¬ 
stroying QCPs exist in TMOs. 

A Simple Perturbed Kitaev Model 

Consider the Kitaev model in an external Zeeman mag¬ 
netic field: 

^ E -^-E^" 4) 

<hj> i 

Using the D = 2 Jordan-Wigner mapping [9], H{see sup¬ 
plementary informations [? ]) can be brought to the fol¬ 
lowing instructive form (this also bares the exaet emer¬ 
gent low-dimensional gauge symmetries, as pointed out 
by Fradkin et al. [10]): H = Hk + Hz, where 

hk + ^( 4*^4 + ^- 4 ] 

Q 

+ ~ ~ 1 ) ( 2 ) 

and Hz = ‘^hz^^{clai + h.c). This is just a spinless 
Falicov-Kimball model with a local hybridization hz in 
d = 2, but with an additional ‘spin’-triplet pairing term, 
whence topologically ordered ground state of Hk now 
appears as a spin-triplet BCS superfluid of JW fermions. 
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Here, we define a'- = (2n^Q, — 1) = iBiiBi 2 = dzl, which is 
now a dynamic gauge field variable on the center of each 
zz-links of the honeycomb Kitaev model (when Ui^cx — 
0,1 for all i. This transforms the original spin model on 
a honeycomb lattice to a two- “orbital” spinless model 
on an effective square lattice. While [nio^^H] = 0 for 
all i (reminiscent of the spinless FK model) ensures the 
exact solvability of the KM when = 0, a finite Zeeman 
field immediately spoils the beautiful integrability of H. 
Interestingly, however, finite hz simply translates to a 
local hybridization between the c, a JW fermions. These 
observations are crucial for an ‘almost exact’ analysis of 
the KM in a Zeeman field as detailed below. 

Alternatively, we motivate the instabilities of the TO 
state of Hk as follows. In Kitaev’s language, a par¬ 
tial local Z 2 symmetry associated with [b^bj^HK] = 0 
for all {ij) and a = x^y still holds when hz 0 (only 
[6f Hk] 7 ^ 0 along zz-links, since Hz = i mixes 

6^, c on each site). This directly implies a dimensional re¬ 
duction and associated lower-T)(= d < D) (here, d = 1) 
gauge-like symmetries (GLS) naturally emerge. The only 
order that survives the action of these GLS is a direc¬ 
tional spin-nema^zc order [10], whose order parameter, 
(n) = ) is the only one that remains 

invariant under the emergent d = 1 GLS: thus, the TO 
of Hk is only partially “melted” by a Zeeman field. The 
ultra-short-ranged spin-liquid phase known for hz = 0 
thus immediately becomes unstable to a deconfined crit¬ 
ical liquid inherently unstable to incipient spin-nematic 
order (on xx,yy zig-zag links) co-existing with a (with fi¬ 
nite fiux on zz-links) field-induced zz-spin order, for any 
hz 7 ^ 0 dit T = 01 It also turns out that the zz-spin 
correlation function, Xzzi'l — j) — {'I — (iiof shown, 
proof essentially identical to Tikhonov et al [11]) using 
a simple bubble composed of the {i ~ j)^ Gcd'l — j) 
propagators for H in Major ana language. Remarkably, 
field-induced magnetization now appears as an excitonic 
condensate of JW-fermions. 

In JW-fermion language, for special values of hz, the 
now dynamical Z 2 flux can, by itself, self-organize into 
stripe-ordered patterns, in analogy with crystal states 
found in the FK model [12] (here, Uia = (1 + (ai)/2 play 
the role of immobile fermions in the FK-like model for 
hz =0). In the present case, such crystals occurring with 
a large range of periodic structures have a novel origin as 
condensation of topological kink dipoles (a condensation 
of {iBiiBi 2 ) on each zz-bond). On the other hand, field- 
induced metamagnetic transitions, now associated with 
‘excitonic solid’ states, are also generically expected from 
a Falicov-Kimball model with hybridization. 

We analyze H within dynamical mean-field theory 
(DMFT). Our choice is motivated by features special 
for Hk (i) when hz = 0, the KM possesses ultra¬ 
short-range spin correlations [13] rigorously uncorre¬ 
lated beyond nearest neighbors, (ii) a DMFT-related ap¬ 
proach [14] thus gives the first instance of an exact dy- 




(d) 


FIG. 1. Upper panel: Magnetization and static susceptibil¬ 
ity as a function of externally applied field for (a) Jz = 0 , 
(b) Jz=0.2Jx, and (c) Jz^OAJx. Lower panel: Plateaus and 
jumps in the anisotropic limit of the KM, with Jz = 0.67^ 
(d). In the inset the 77 is shown as a function of {Jx = 0.5). 
The axes labels for figures (a), (b) and (c) are same as (d). 

namical spin fluctuation spectrum in a I) = 2 model, and 
(hi) the fact that spin correlations on all length scales 
along XX, yy bonds have been exactly subsumed into one- 
fermion hopping and p-wave BGS pairing terms obviates 
the usual difficulty of ignoring long-range spatial correla¬ 
tions, wherein any mean-field theory, including DMFT, 
is usually critiqued in D = 2. Moreover, Hz also has a 
local excitonic form and the zz-couplings appear as lo¬ 
cal interaction. Thus, these novelties, again specific for 
the KM, make the DMFT almost exact and allows for 
exhaustive analysis. 

The impurity model of DMFT is solved by multi¬ 
orbital iterated perturbation theory as an impurity 
solver: though not “exact”, its proven quantitative 
accuracy vis-a-vis exact QMG results for the spinful 
Anderson lattice model [15] ensures its efficacy in this 
case. We focus on the c, a-fermion spectral functions 
and field-induced magnetization, now simply equal to 
a local excitonic average m{hz) = {{c[ai h.c)). To 
facilitate discussion, we show the evolution of m{hz) 
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FIG. 2. Plateaus and jumps in the ‘magnetization’ in the 
isotropic limit of the KM. The corresponding susceptibilities 
are also shown. 


vs. hz in Fig. 1. When Jz « Jx = Jyi i^ihz) rises 
monotonically with hz, as expected of a field-polarized 
paramagnetic phase. However, its underlying topological 
origin is underscored by observing that m(hz) — h'2 
with 1.0 > ?7 > 0.85 for 0 < < 0.25(= Jxl^)'^ for 

larger Jz, r]{Jz) reverses its slope (Fig. 1 lower panel 
(inset)), increasing monotonically toward 1.0. The 
special character of JzjJx — 1/4 is related to its being 
precisely twice the (mean-field) excitation energy of 
the three fiat bands that arise in a mean-field theory 
(related to the original resonating valence bond (RVB) 
description in the cuprate context [16]) distinct to 
ours [17]. It marks the point where the a-fermion level 
and concomitant excitonic tendency begins to play a 
role for JzjJx ^ 1/4. That m(hz) starts developing 
plateaus (or incompressible JW-excitonic solid phases) 
as a function of hz only for Jzj Jx > 1/4 also supports 
this view. 

Remarkably, for JzjJx > 1/4, m{hz) as a function 
of hz shows an extensive sequence of magnetization 
plateaus, now associated with incompressible excitonic 
solid phases, separated by excitonic “fluid” phases. 
These are intimately linked to the kink-dipole “con¬ 
densates” (crystals of the “flux”, {iBiiBi 2 ): see the 
'o^a{hz) VS hz in Suppl. Info). These are thus first- 
order quantum phase transitions between sequences of 
“topological crystal” states, each sandwiched between 
topological Bose-Einstein condensates of excitonic kink 
dipoles. These are nothing but zeeman-field induced 
metamagnetic transitions with accompanying hysteresis 
at each magnetization jump (see Fig. (3b) of supplemen¬ 


tary informations), and go hand-in-hand with jumps in 
{i^a{hz) at the same hz values. These topological crystal 
phases occur (Fig. 2) for a range of odd-denominator 
valued magnetizations, as well as at certain even de¬ 
nominators (relative to saturation, at m{hz)/rusat = 
1/16,1/12,1/11,1/9,1/8,1/7,2/13,1/6, 2/11,1/5, 2/9, 
1/4, 2/7,1/3, 3/8, 3/7, 2/5,4/7, 5/8, 3/4,...). Odd denom¬ 
inator plateaus in the off-diagonal conductivity are well- 
known in the fractional quantum Hall effect (FQHE), 
which is also the only other example of a real system 
rigorously exhibiting topological order. In spin systems, 
such exotica have also been found and widely discussed 
in frustrated Shastry-Sutherland models [18, 19], where 
they emerge due to competition between frustration and 
field-induced magnetism (even-denominator plateaus 
are also visible there). So our plateaus are not related 
to EQHE-like physics, even though (in contrast to 
earlier work) they arise in a perturbed (field-induced) 
partially TO phase. Nevertheless, in Eig. 2, we exhibit 
m{hz) and xi^z) = {d/dhz)m{hz) together to illustrate 
“analogies” with EQHE. We interpret the oscillations 
in x{dz) ns “de Haas van-Alphen (dHvA)” quantum 
oscillations of topological JW fermions in a partially 
magnetized spin liquid phase. Such exotica may show 
up [20] in magnetic torque studies. This reflects “hid¬ 
den” coherence in a spin liquid: interestingly, hidden 
off-diagonal long-range order (ODLRO) was conjectured 
by early on by Anderson precisely in this context [21]. 
Here, hz 0 partially polarizes the spinons reflecting in 
magnetization plateaus, and these ordered phases “ ride 
on top” of a partially topological ordered state. This 
component has in-built coherence, thanks to the p-wave 
BCS term in manifesting as coherent “quantum 

oscillations” we unearth. Thus, in stark contrast to 
conventional paramagnets, the TO phase of the Kitaev 
model undergoes an intricate sequence of partial ordering 
“solid” instabilities, co-existing with the remnant of the 
TO fractionalized spin liquid in an external Zeeman field 
before reaching full saturation. 

The DMET spectral functions reveal deeper insights 
(See supplementary informations). Clear “orbital” se¬ 
lectivity is seen in c and a-fermion spectral functions, 
and this is a direct fall-out of the in-built differentiation 
between these states in the JW-fermionized Hk itself. 
While hybridization and local Hubbard correlations con¬ 
spire to produce a JW “Kondo” insulator for a-fermions, 
they produce large-scale spectral changes in the c-fermion 
sector. Eocussing on the variations in pc{oo), Pa{^) across 
representative samples where magnetization steps ap¬ 
pear, it is clear that large-scale reshuffling of dynamical 
spectral weight across energy scales of order the band¬ 
width (equal to 2 Jx here) accompany emergence of each 
plateau 

Thus, the plateau structure has its origin in the ubiq¬ 
uitous competition between Mott localization (Hubbard 
term) and itinerance (hybridization and hopping) in iJ, 
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(a) 



FIG. 3. The renormalized dispersions ((a) and (b)), corre¬ 
sponding Luttinger surfaces ((c) and (d)) for dispersive JW 
Fermions at plateaus 1/3 and 2/3 respectively. The bound¬ 
ary of yellow and blue regions corresponds to the contour of 
ky) = 0 and has been marked similarly inside the figures 
in the lower panel, kx and ky are in units of 7r/a, where a in 
lattice constant. 


and the field simply tunes the balance between these ten¬ 
dencies, now by favoring the tendency to exciton forma¬ 
tion. This is reminiscent of the state of affairs character¬ 
izing the FL* view proposed in the context of /-electron 
quantum criticality [22]. 

More unanticipated features show up in DMFT “band 
structures”. In Fig. SI, we exhibit the c-fermion dis¬ 
persions and “Fermi” surfaces obtained as usual from 
the DMFT propagators. Remarkably, we find that 
the “Fermi surface” corresponds to zeros^ rather than 


poles, of Gcc{k,uj): this is evident from the structure of 
Gcc{k^uj) itself, which features an anomalous “BCS”-like 

self-energy structure given by T.an{k,uj) = 
so that the poles of T^an appear as zeros of Gcc{k^uj). 
Thus, this is a generalized Luttinger surface, rather 
than a conventional Fermi surface. This is analogous to 
what happens in correlation-driven selective-Mott tran¬ 
sitions in /-electron systems in the context of the FL* 
phase, or in the YRZ ansatz [23]. In our case, clear 
changes in topology of the Luttinger surface occur at each 
plateau, and are thus intimately linked to correlation- 
driven Lifshitz-like transitions accompanying each field- 
induced plateau. Moreover, given that the c-fermion 
spectral function is that of p-wave nodal Bogoliubov 
quasiparticles, this this state is a p-wave nodal spin semi¬ 
metal. These novelties are signatures of the underlying 
partial TO surviving in the partially magnetized phases. 
A Zeeman field results in the transfer of dynamical spec¬ 
tral weight from the fractionalized spin liquid component 
to the excitonic solid, but it is indeed remarkable, thanks 
to the undisturbed emergent d = 1-GLS, that it partially 
preserves the original TO. 

Finally, the large-limit of Hk is known to be noth¬ 
ing but the Kitaev Toric Code (TC) model [1]. Here, this 
implies a Mott-\\ke insulating phase of the JW fermions. 
The no double-occupancy constraint of a c, a fermion 
on each site i is naturally implemented by a Gutzwiller 
projector, Pq = n^(l — nicUi^) acting on the p-wave 
BCS ground state. Explicitly, |4 ^)tc = Pc^ki^k + 
'^kclc^_j.)\vac) where \vac) is the c-fermion vacuum. But 
this is precisely the Gutzwiller-projected (p-wave) BCS 
state, long known as one of the optimal states describing 
a (d-wave) RVB spin liquid state in the context of the 
famed RVB description of the cuprates [16]. Thus, the 
ground state of the TC model is just the p-wave RVB 
state! This suggests the enticing possibility of investigat¬ 
ing TC physics and its non-Abelian excitations in future 
in terms hitherto developed for quantum spin liquids in 
quantum Heisenberg models [24]. 

We now detail a specific quasi-realistic scenario where 
our findings could aid finding of novel orders in a real- 
world set-up: 

1. Kitaev Josephson Junction Arrays As men¬ 
tioned before (see also supplementary informations), a 
novel implementation of the perturbed KM using JJA 
has been proposed [3]. Our results can now be used to 
unearth novelties in such engineered systems. 

We now see that varying the flux through a qubit 
leads to a generation of a whole host of ordered “JW- 
excitonic cooper-pair” crystals, with a succession of peri¬ 
ods set by the competition between hz (flux) and the ma¬ 
terial parameters of each JJ (these enter the Kitaev cou¬ 
plings [3] in Jx^Jy = These are now interpreted 

as kink-dipole crystals, and arise as an explicit mani- 
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festation of the fact that strong geometric frustration- 
induced topological liquid state of the KM undergoes an 
intricate sequence of “partially ordered” solid phases be¬ 
fore achieving a full solid at saturation. In view of the 
co-existing and remnant topological critical spin-liquid 
arising from emergent “d = 1” GLS in this state, one 
may also refer to each of these states as critical topo¬ 
logical supersolids: these turn out to be unstable to ne¬ 
matic order ^Jx^Jy (see supplementary informations). 
A direct manifestation of this insight is that, in anal¬ 
ogy with excitonic Josephson effects in electron-hole bi¬ 
layers [25], one expects the Josephson critical current, 
Jc, to depend on the excitonic correlations. Explicitly, 
Jc ^ hKclai] c^^aj) hKclai) {c^^aj). As a result of the 
p-wave “BCS” term in the fermionized Hx, there will 
also be a direct “fermionic” current, now caused by the 
p-wave JW-fermions. Thus, remarkably, we expect the 
critical “Josephson” current to exhibit fractional oscilla¬ 
tions as the flux is slowly ramped up. Fractional Shapiro 
steps are also directly predicted as a consequence of the 
underlying superstructure of kink-dipoles. Inducing tran¬ 
sitions between successive steps with a pulsed electro¬ 
magnetic field can result in emission of coherent electro¬ 
magnetic radiation whose energy hft can be tuned in this 
setup. With realistic JJ parameters, ft can be tuned to 
he in the GHz or sub-GHz region: thus, this setup can 
achieve flux-tuned control of EM radiation. Excitations 
of such ordered states as those found here will be slid¬ 
ing collective modes of the excitonic solid akin to phason 
modes in the GDW context, and these are precisely the 
“Josephson plasmons”, suggesting novel applications to 
“topological plasmonics”. Ours is thus a first proposal to 
exploit an almost exactly solvable model of fractionalized 
JW-fermions in = 2 to detect topological analogues of 
famed quantum interference effects. Other subtle man¬ 
ifestations of ODLRO, such as topological analogues of 
Little-Parks experiment, etc, could also be imagined in 
the JJA context. Further, to the extent that it is difficult 
to generate single JW fermions by splitting up JW exci- 
tons, the states we find have a much more efficient topo¬ 
logical protection (the topological liquid phases are in¬ 
herently symmetry-protected topological (SPT) phases), 
these findings can be exploited for possible topological 
quantum computation (TQG) engineering, along lines 
proposed by You et al [3] . 

Conclusion 

To conclude, we have presented strong evidence link¬ 
ing a large sequence of “excitonic ordered solid” phases 
of the Zeeman-held driven KM to novel instabilities of a 
rigorously known TO phase. These could be touted as 
topological analogues of the famous ”Barkhausen steps” 
encountered in the magnetization process of a conven¬ 
tional held-driven magnet. Here, however, they corre¬ 
spond to partially magnetically ordered complex unit¬ 
cell patterns emerging as partial ordered states of a frac¬ 
tionalized spin liquid, betraying their highly unconven¬ 


tional nature. Suitably engineered JJAs as proposed [3] 
could unearth the structures we propose by hux tun¬ 
ing, with novel applications to novel quantum interfer¬ 
ence phenomena, plasmonics and TQG. Finally, if ex¬ 
amples of such perturbed KMs to real transition-metal 
oxides could be found, manifestations of the underlying 
TO emerging as lower-d(= 1) gauge-like symmetries, to¬ 
gether with their intimate link to unconventional nematic 
order, could reveal themselves under appropriate exter¬ 
nal perturbations. This remains an enticing perspective 
for future work. 
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Supplementary Informations 

D = 2 Duality Transformations 

Kitaev’s famed solution relies on mapping spin-1/2s to 
bilinear s of four ^ b'^ ^ ^ c) Major ana fermions along 
the three (xx^yy^zz) bonds (Fig. SI). We start with 
an alternative but equivalent approach due to Chen and 
Hu [S 2 ] and Nussinov [S3] that exploits a duality map¬ 
ping of the spins via a D = 2 Jordan-Wigner (JW) 
fermionization of H. Remarkably, since this affords an 
exact mapping between topological and classical orders 
in D = 2 spin systems, studying TO(s), their excitations 
and instabilities is transmuted into the better-known lan¬ 
guage of well-known models hosting classical order(s) 
characterized by symmetry-breaking. Thus, 

(SI) 

and 

St^ = {2c\^Cij - 1) (S2) 

and defining Majorana fermions on the two (“black” and 
“white”) sublattice sites of the honeycomb or brick-wall 
lattice as = (c — = (c -h c^)yj and = 

(c + = (c — c^)b/h followed by the introduction 

of fermions c= {Ayj iAb)/2^c^ = {Ayj —iAb)l2. Using 
these transformations, we get the fermionized form of 
H = used in the main text. In H as thus defined 

in the main text, a[ = iBi\Bi 2 is a bilinear of Majorana 
fermions. Since = 0 for all i, it takes eigenvalues 

±1. Thus, we define a'- = (2n^^Q, “ 1)^ where rii^a = O 5 1 
and [rii^a^ H] = 0 for all i as well, preserving the local Z 2 
symmetry . Using this representation, the zeeman term 
maps onto the excitonic term in JW-fermion language. 

Non-singular Spin Fluctuation Spectrum 

Here, we point out how the “Dirac” spectrum for c- 
fermions cuts-off the infra-red (X-ray edge) singularities 
in the Kitaev model with hz = 0. In a normal Fermi 
liquid, the infra-red singular nature of the localized prop¬ 
agator (as well as of the “excitonic” propagator) can be 
traced back to the appearance of logarithmic divergences 
in the particle-hole polarization function. This necessi¬ 
tates summing of infinite number of Feynman diagrams 
for any arbitrary coupling strength [S4], or “summing 
the logs”, leading to a power-law singularity in dynam¬ 
ical response functions. In the Kitaev model, however, 
two features mitigate against this: first, the immobile- 
fermion level lies at —Jzl‘ 2 . rather than at u; = 0 , and, 
second, the “Dirac”-like DOS for the c-fermions (thanks 
to the p-wave BCS term in Hk) cuts off the infra-red sin- 
gularities. Explicitly, writing G'c«a|(o)(ic^) = iu,+j,/2 
Gcc{k,iUj) = 

near the ‘Dirac” points, we find that the one-loop po¬ 
larization function. 


nO)r ] = - f ^ f _ 1 

J 271 J (27iY {ii/Y — {vkY iv — iuj ^ Jz/2 

(S3) 

behaves as n^^^(ici;) (27ru^) ^ Jz/ 2 )hi{{Jz/ 2 — 

iuj)/W) — VF], where IF ^ is a cut-off of order the 
spinon bandwidth. The log-singularity is now cut-off by 
the factor {iuj ^ Jz/ 2 \ which itself arises from the linear 
DOS of the non-interacting c-fermions. Thus, there is no 
infra-red singularity in this case, and, as expected from 
this line of argument, a RPA calculation using the bare 
propagators indeed captures the essential features of the 
exact solution to a high accuracy [S5]. 

In the TO phase {hz = 0), the c-fermions thus expe¬ 
rience a strongly fluctuating {Z 2 ) potential, which ‘sud¬ 
denly’ switches between 0, Jz during every hop of the 
c-fermions. If the pairing term were absent, this would 
constitute a complete reconstruction of the entire Fermi 
sea due to the Anderson-Nozieres-de-Dominicis “orthogo¬ 
nality catastrophe” (OC). Along with the rigorous ultra- 
short spin correlations in Hk^ this is what makes the 
‘local’ approximation exact here. But with hz = 0, Cq 
has a Dirac-like form which cuts off the infra-red singu¬ 
lar behavior and the OC, since the DOS vanishes linearly 
at cj = 0(= ei?) as shown above; this will not be true any¬ 
more when hz is finite and Fermi pockets, corresponding 
to finite DOS at cj = 0, are expected to appear in the 
c-fermion spectrum. However, this is still not enough, as 
the local hybridization also produces finite recoil of the 
(a-fermion: this must cut-off the X-ray-edge singularities 
in spectral functions. But since any finite hz also leads 
to an excitonic condensate at T = 0, this provides an 
additional route to removal of infra-red singularities by 
a (partial in our case) coherence-restoring instability to 
field-induced order. This is precisely what we find (see 
main text). 


Non-Trivial Topological Order with hz ^0 

Here, we show how the emergent D = 1 GLS men¬ 
tioned in the text are linked to a novel “hidden” topo¬ 
logical order even when hz 0. 

Noticing that a Zeeman field along ^^-bonds implies 
[ib^b^^Hx + Hz] = 0 for all {ij)\\xx^yy implies persis¬ 
tence of partial local Z 2 gauge symmetry and emergent 
low {d = 1) gauge-like symmetries, we show that these 
hide a novel kind of topological quantum critical point 
(TQCP), separating two topological phases. Since topo¬ 
logical order can be described by conventional (Landau¬ 
like) order when expressed in terms of dual (Jordan- 
Wigner) variables, we take this approach. Remarkably, 
it will turn out that the special point Jx = Jy{= J) sep¬ 
arates two IsingAdke nematic phases at T = 0. 

In the single chain case, Hk is trivially diagonal¬ 
ized to give two quasiparticle bands with dispersions 
J- + 2 JxJyCos{k) (with one band fully 

occupied and the other empty). The second derivative 
of the ground state energy diverges at Jx = Jy, where 
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FIG. SI. (Color online) JJA on a honeycomb lattice( [SI]). 


the excitation gap vanishes. This quantum transition at 
T = 0 does not involve any symmetry change, rather, it 
is a change in topological order. This is made explicit by 
the transformation [S6, S7] 

= rf_irf,5f (S4) 

which recasts Hk as the D = 1 quantum Ising model 
(QIM) on the dual lattice. When > J^the dual spin- 
correlation function shows long-range order: 

Limi^oo{ToT2i) ~ [1 - (S5) 

In the original spin-representation, this is a string order 
with “order” parameter Ax{i) = 

For Jy > Ja,, we exploit the well-known Kramers- 
Wannier duality of the D = 1 QIM to conclude that the 
order is that of strings of the 5'^, i.e, the order parameter 
is A, = 

At Jx = Jy, both these order parameters collapse to 
zero simultaneously. The topological excitations which 
destroy each of these orders are precisely proliferated 
string-like excitations associated with the dual order: 
these are consequences of the emergent d = 1 GLS. Re¬ 
markably, it is easily seen that the two non-zero order 
parameters are nothing but the only ones invariant un¬ 
der the emergent GLS, i.e, they are the spin-nematic or¬ 
der parameters, defined by (n) = {SfSj — Sf S^). Thus, 
Jx = Jy is a TQCP separating two topological ordered 
phases. In dual language, it separates two Ising-nematic 
ordered phases. 


Implementation of the Perturbed Kitaev Model 
in Engineered Josephson Junction Arrays 

Here, we briefly recapitulate the proposal for simulat¬ 
ing the perturbed Kitaev model we have studied by us¬ 
ing a specific engineered Josephson junction array on the 
honeycomb lattice. 

Our results are directly applicable to the engineered 
JJA of You et al. [SI] as long as we identify the Kitaev 
spins as pseudospin-1/2s associated with a pair charge 
qubit: defining | | as two charge states with cooper 

pair number Ui = 1,0, the S' = 1/2 representation per¬ 
mits a rigorous mapping, m = {1 — Sf)l2^ cos0^ = Sf/2 
and sin^i = —S^ 12 with (j)i being the pair phase at i, 
conjugate to the number as implied by Ui = —If 
each charge qubit is placed at an “optimal” point where 
Ug = 1/2, and Jx{= Jy) < Jz is considered, one ends 
up with the anisotropic Kitaev model in an “external 
field” h^-Sf {not hSf). However, this is of no con¬ 
sequence, as we obtain our model simply by relabelling 
Sf ^ —Sf' Further, we then have J^ < Jx{= 

J), which is precisely the limit we are concerned with. 
Remarkably, since h = Eji{^i) = 2J^jcos(7r^i/^o) with 
^0 = 7rfi/e the flux quantum in the JJA, changing the 
Josephson coupling energy offers a subtle knob to tune 
the “Zeeman” field in the KM, and, in fact, in the JJA 
context, a real external magnetic field or adding non¬ 
magnetic impurities can change Eji{^) [S8]. 




FIG. S2. (Golor online) The DOS of the interacting c and a 
(left and right respectively) fermions at plateaus | and | . 


DMFT Spectral Functions 

Here, we expound briefly upon our DMFT results for 
the spectral functions used to discuss spectral features 
and their consequences in the main text. 

We have used the two-orbital iterated perturbation 
theory (IPT) as an impurity solver in the context of solv¬ 
ing H using DMFT. Though not “exact” in the sense 
of continuous-time quantum Monte Carlo (CT-QMC) 
or the numerical renormalization group (NRG) impurity 
solvers, its numerical efficacy in the context of DMFT 
for a related spinful Anderson lattice problem vis-a-vis 
CTQMC is well-documented, allowing a good degree of 
confidence for its use in our case. 
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FIG. S3. (Color online)(a) na{hz)lna,sat as a function of 
magnetic field (hz) in the isotropic limit (ria^sat is the sat¬ 
uration value of ria). The jumps and plateaus in the profile 
closely correspond to the same numbers in nn(hz)Irrisat in the 
same limit, (b) Hysteresis in magnetization. 

In Fig. S2, we show the c- and the a-JW fermion 
spectral functions for the isotropic Kitaev case at values 
m{hz)lmsat = 1/3, 2/3. Several notable features stand 
out clearly. First, the in-built “orbital” differentiation in¬ 
volved in the pure KM has direct consequences on their 
respective evolution as a function of hz. In particular, 
the (a-fermion DOS resembles that of a Kondo insulator, 
owing to the fact that a combination of hz induced c — a 
hybridization and the bare energy level of the a-fermion 
in the presence of a “Hubbard” (Jz) interaction term 
conspire to generate such a state for the parameter val¬ 
ues set by H. However, the c-fermion DOS shows drastic 
spectral changes as a direct consequence of the dynami¬ 
cal interplay between the “Dirac” band structure of the 
c-fermions, c — a hybridization, and the Hubbard interac¬ 
tion. Notwithstanding this, however, the c-JW fermions 
continue to be elementary excitations of the perturbed 
spin liquid, since ImScc(co’) always vanishes at low en¬ 


ergy in DMFT (not shown). This aspect is a direct conse¬ 
quence of the important observation, namely, that a finite 
hz leads to emergent lower-dimensional d{= 1) < D = 2 
gauge-like symmetries: thus, these spinon-like excitations 
are topological excitations (non-local strings in the orig¬ 
inal spin variables in Hk) of this partially TO state. 
However, this fractionalized component co-exists with a 
finite “excitonic” order parameter, which is just the field- 
induced magnetization (Fig. S3 (a)) in original spin lan¬ 
guage. 

This has important consequences: clearly, our results 
show that, across successive magnetization plateaus, 
large reshuffling of dynamical spectral weight over scales 
of order the full c-fermion band-width occur. This clearly 
attests to the important role of dynamical correlations in 
this interplay and generation of the plateau structures. 
To the extent such dynamical correlations are crucial, 
we opine that less sophisticated approximations such as 
static Hartree-Fock may miss important aspects of the 
detailed structures found here. 
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